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The goal of this paper is to highlight several issues which are most crucial for the understanding 
of the "metal-insulator transition" in two dimensions. We discuss some common problems in in- 
terpreting experimental results on high mobility Si MOSFETs. We analyze concepts and methods 
used to determine the critical density of electrons at the metal-insulator transition. In particular, we 
discuss the origin of the temperature dependence of the resistivity and reasons for this dependence to 
flatten out at some electron density in the vicinity of the metal-insulator transition. This flattening 
has recently been proposed to indicate a true quantum phase transition. We suggest an alternative 
interpretation of this result and demonstrate the consistency of our proposition with the experimen- 
tal data. One of the main questions, which arise in connection with the transition, is whether or not 
the metallic state is qualitatively distinct from a conventional disordered Fermi liquid. We analyze 
the arguments in favor of both affirmative and negative answers to this question and conclude that 
the experimental results accumulated up-to-date do not provide convincing evidence for the new 
state of matter characterized by a metallic-like residual conductivity. We also discuss in details the 
measurement and control of the electron temperature; these issues are crucial for interpreting the 
low-temperature experimental data. 



PACS numbers: 71.30.+h, 72.15 Rn, 73.40. Qv 
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I. INTRODUCTION 

After several years of intensive experimental and the- 
oretical efforts (see, e.g., Ref. [ [l] for an extensive bib- 
liography), even the basic features of the phenomenon 
known as the "metal-insulator transition in two dimen- 
sions" (2D MIT) remain to be the subjects of ongoing 
discussions. Is this phenomenon a true "quantum phase 
transition" (QPT) [|] or can it be understood in terms 
of conventional physics of disordered conductors [ MM ? 
This question is at the heart of the whole discussion. 
Recently, we wrote a paper [ [| in favor of the second 
possibility. In particular, we argued that 

i) in the "metallic phase" , 2D systems seem to behave 
as a quite conventional disordered metal rather 
than a distinctly new state of matter and 

ii) it is possible to explain the anomalous behavior of 
the resistivity, p(T), by an interplay of two temper- 
ature dependences: the one given by a metallic-like, 
quasiclassical (Drude) resistivity [ ||, Pd{T), and 
the other one arising from quantum localization. 

We do realize that ii) represents a rather naive ap- 
proach, at least because it does not fully take into account 
the Coulomb interaction between electrons, whereas the 
most pronounced effects have been observed in systems 
with a priori strong electron-electron correlations. Nev- 
ertheless, this approach allows one to describe qualita- 
tively, and even semi-quantitatively, most of the results 
on electron transport both in metallic phase and near the 
transition point. This suggests that it is a good starting 
point for developing an adequate theoretical understand- 
ing of the MIT in two dimensions. 

Conclusion i) is based on the analysis of existing ex- 
perimental data, as summarized in Refs. [ ^,0. It has 
been further supported by recent experimental papers [ 
p],p|"|Tot, which showed that the transport properties of 
Si/Ge, Si MOSFET and p-GaAs structures in the metal- 
lic regime can be successfully interpreted in conventional 
terms. Authors of Refs. [ ^,|l0| have also identified the 
contribution of electron-electron interactions to the resis- 
tance (via measuring the temperature-dependent part of 
the Hall resistance) and found that it remains quite small 
even when parameter r s is rather large. There is also 
a number of recent theoretical papers [ prj-pTP] which, 
although differing in a particular mechanism for the In- 
dependence of the Drude resistivity, share the general 
spirit of propositions i) and ii). 

An alternative point of view, expressed, for example, 
in Refs. [ fl5-17|, is that the observed MIT-like behavior 
indicates a true quantum phase transition between an 
insulator and a novel metallic phase. Motivated by the 
importance of this controversy, i.e., QPT versus proposi- 
tions i) and ii) , for the field of quantum transport in 2D 
and its possible relevance for other realizations of quan- 
tum phase transitions, we decided to analyze in detail the 
arguments on both sides of the QPT question. In this 



paper, we discuss recent and some of the previously pub- 
lished experimental results on the resistivity of Si MOS- 
FETs in the vicinity of the MIT, along with some of the 
theoretical interpretations of these results. We conclude 
that that there is no convincing experimental evidence 
for the QPT in the existing data. 

As one of the main questions in the field is the behavior 
of a 2D system in the limit T — > 0, we begin our analysis 
with a short summary of recent results by Kravchcnko 
and Klapwijk (KK) [[jj, who measured the temperature 
dependence of the resistivity down to a bath temperature 
as low as Tbath = 35 mK. This paper reports measure- 
ments on a single Si MOSFET sample with peak mobil- 
ity /ipeak = 27,000cm 2 /Vs in the range from 35 mK to 
1.2 K at five different electron concentrations. Digitized 
data from Ref. [ |T^] is shown in Fig. [l] (curves 1-5). 

In Ref. [ [l5|] , the following two points are emphasized: 

1) For n = n\ = 6.85 x 10 10 cm~ 2 and n = n,2 = 
7.17 x 10 10 cm~ 2 (curves 1 and 2), the resistivity 
decreases as temperature increases {dp/dT < 0); 
for n = Hi — 7.57 x 10 10 cm -2 and n = n 5 — 7.85 x 
10 10 cm -2 (curves 4 and 5), the resistivity increases 
with temperature (dp/dT > 0). Based on these ob- 
servations, the authors argue that curves 1 and 2 
correspond to an insulating phase, whereas curves 
4 and 5 demonstrate a metallic behavior. At n = 
n 3 = 7.25 x 10 10 cur 2 (curve 3), only a weak (±5%) 
temperature dependence was observed within the 
interval of bath temperature from 35 mK to 1 K. 
The authors conclude that nz = 7.25 x 10 10 cm -2 is 
the critical electron density for their sample, i.e., it 
corresponds to the metal-insulator transition point. 
Note that this conclusion is based entirely on the 
temperature independence of curve 3 over a rela- 
tively narrow temperature range. 

2) Another observation, emphasized in Ref. [ |T^], is 
that in the range 0.1 K< T <0.4K the metal- 
lic p(T) - dependences (curves 4> 5) appear to be 
nearly linear and thus different from the exponen- 
tial temperature dependence observed in other ex- 
periments. KK stress that in this region the re- 
sistivity shows no insulating up-turn at the lowest 
temperature achieved. 

KK consider the fact that they observed no tempera- 
ture dependence of the resistivity of their sample at a cer- 
tain concentration n = 713 as a strong argument against 
our model [ ^ . They point out that such a precise cancel- 
lation would require a special relation between Pd(T) and 
the scaling /3(p)-function. KK argue that this relation be- 
tween two objects of entirely different origin "would be 
a remarkable coincidence" , and is therefore unlikely. 

Based on the arguments listed above, authors of Ref. 
[ |l5| conclude that their experiments are consistent with 
the existence of the zero-temperature quantum phase 
transition . 
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FIG. 1. Resistivity vs temperature, reproduced 
from Fig. 2 of Ref. [15]. Electron density n in 
units of 10 10 cm~ 2 : 1 - 6.85, 2 - 7.17, 3 - 7.25, 
4 - 7.57, 5 - 7.85. Dashed curves Ai and A^: 
fits of data 1 and by 2.29exp [(403/T) 1/2 ] and 
2.84exp [(11.11/T) 1/2 ], respectively (T is in mK). 
Solid curves B\ and B2: fits of the same data by 
1.25 exp [(1685/T) 1/3 ] and 2.55 exp [(8./T) 1/3 ] , re- 
spectively. Semi-open symbols (1"): reflection of data 
points 1 about curve 3. Curves E and F: best fits of 
data 4 and 5 by p A [T) = 1.2 + 1.45 exp(-450/T) and 
p 5 (T) = 0.63 + 1.6exp(-485/T), respectively. 

Note that a rather narrow range of densities was ex- 
plored in Ref. [ [l^] . To show the development of the MIT- 
phenomenon over a wider range of densities and temper- 
atures, we present in Fig. || the results of Ref. [ |l8| for a 
sample with a close value of /i pca k- The encircled region 
of the {T, 71} plane (region KK in this figure) indicates 
roughly the domain of parameters explored in Ref. [ [l5| . 
The temperature in Fig. 2 is normalized to the Fermi en- 
ergy Ep , in order to demonstrate that resistivities corre- 
sponding to quite different Fermi energies, and therefore 
r s -values, exhibit similar T-dcpcndences. A strong (50- 
100%) drop in p is still present for n < 20 x 10 n cm~ 2 , 
which corresponds to r s ~ 2 and p ~ 0.01/i/e 2 . Thus, 



contrary to the popular opinion, the resistance drop is not 
intrinsic only to the low density range (high revalues). 
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FIG. 2. Resistivity vs temperature (from 0.3 to 
~ 40 K) for a wide range of densities and Fermi 
energies. The figure is reproduced from Ref. [18]. 
Encircled region: domain of parameters explored in 
Ref. [15]. Dash-dotted vertical line depicts the em- 
pirical temperature, Tq — 0.007-Ef, below which the 
logarithmic temperature dependence sets in. 

We would like to stress once again that in the region 
of high densities and small resistivities quantum inter- 
ference corrections of any kind (weak localization, ex- 
change and Hartree interactions, spin-orbit, etc.) would 
amount only to a few percent variation (of both signs) 
in p with T. The phonon contribution to the resistiv- 
ity of Si-MOSFETs can be shown to be negligible at 
least for n < 20 x 10 11 cm~ 2 in the relevant T-range. 
Therefore, there must be some other mechanism of the 
T-dependence at work, neither of quantum interference 
nor of phonon nature. The resistance drop of non-phonon 
nature persists up to the highest temperature of these 
measurements, which is ~ 35 — 40 K for each curve in 
Fig. 2 (for conversion of densities into the Fermi-energies, 
use the formula: E F [K}= 7.31 x n[10 u cm" 2 ]). 
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II. METAL, INSULATOR, AND CRITICAL 
POINT. GENERAL DISCUSSION 

Contrary to a common belief, a metallic phase and 
thus a metal-insulator transition can in principle occur 
in two dimensions. For example, in a hypothetical situ- 
ation of non-interacting electrons spin-orbit coupling re- 
sults in resistivity vanishing with temperature, provided 
that disorder is sufficiently weak ( "weak antilocalization" 
[ |||I|,[l9|] ) . However, if disorder is strong enough, it leads 
to localization. The metallic state is stable with respect 
to weak and short-ranged electron-electron interactions. 
The Coulomb interaction between electrons for small r s 
overwhelms weak antilocalization and thus destroys the 
metallic state with zero residual resistivity. Nevertheless, 
the Coulomb interaction itself brings in an antilocalizing 
contribution to the resistivity (triplet channel) . For small 
r s , this contribution is smaller than the one from the sin- 
glet channel, and thus the net effect of interactions is the 
enhancement of localization. It may be the case that the 
triplet channel wins for larger r s and a metallic state be- 
comes possible. (In fact, this scenario follows from the 
RG treatment [ ^p|-^2|.) However, if such a state exists, it 
must be distinct from the conventional "disordered Fermi 
liquid". 

So far, the experimental evidence is in contrast with 
the existence of a distinctly new state (see Refs.[ 
In fact, at higher densities the resistivity exhibits an in- 
sulating up-turn in agreement with the weak localization 
theory. We believe that in the existing experiments the 
true low-temperature asymptotic behavior has not yet 
been achieved for lower densities (in the vicinity of the 
transition). Therefore, one cannot interpret the crossover 
between dp/dT > and dp/dT < at finite tempera- 
tures as a proof for the existence of two distinct phases. 
Nevertheless, it is instructive to adopt such an interpre- 
tation for a moment and to compare the experimentally 
observed p(T) with the commonly accepted picture of 
a quantum phase transition. Let us first recall what is 
known about the transition in 3D. 

We start our analysis by formulating the definitions of 
metal and insulator. Arguing about definitions is not a 
too rational thing to do - a definition cannot be right 
or wrong. We simply point out that for the well-studied 
case of MIT in 3D, the commonly held definitions of both 
phases differ from those adopted by KK. In 3D systems, 
Pinsuiatoijfr > 0) ► oo, whereas p mc tai(r -> 0) is finite 
(see [ |3p|,|22|] and references therein). In fact, dp/dT is 
negative in both the metallic and insulating phases, pro- 
vided that the system is close to the transition point and 
the temperature is low enough. 

For non-interacting particles, the conventional sce- 
nario of the 3D MIT is well supported by the pertur- 
bation theory and renormalization group arguments [ 
|,|,||J22|. According to this scenario, in the metallic 
phase close to the transition p(T) increases as T~ a as 
temperature decreases. The increase of p{T) saturates 
at T < T sat . Upon approaching the critical point, T sat 



tends to zero. Therefore, exactly at the critical density 
the p(T)-dependence is a power law: p{T) ~ T~ a . In 
the insulating phase, the resistivity behaves as exp(T /T) 
(nearest neighbor hopping) or as exp[(To/T) /3 ] (variable 
range hopping) . The exponent a is determined by mech- 
anisms of dephasing. As a result, it is, strictly speak- 
ing, non-universal. Let the dephasing rate scale with 
temperature as 1/t v oc T p . Then assuming that (a) 
one-parameter scaling holds, (b) this parameter is the 
conductance, and (c) p = 1, one arrives at a = 1/3. 
Note that the electron-electron interactions can change 
the value of a. 

Numerous experiments on 3D doped semiconductors 
confirmed the power-law behavior of p in the critical re- 
gion. However, there is still no consensus regarding the 
value of a: both a = 1/2 [ ||-|2|] and a = 1/3 [ |7||28| 
have been reported. It is also possible that a depends 
on whether a semiconductor is compensated or not. (For 
the review of an MIT in 3D see Ref. [ §§.) 

Returning to the MIT in 2D, we note that the assump- 
tion of a temperature-independent resistivity at the crit- 
ical density is as doubtful as the statement that a = 1/3 
in 3D. Indeed, both of these points can be justified only 
within the one- parameter scaling. This scaling does not 
seem to apply universally even in the 3D case. For a 2D 
system of noninteracting electrons, this scaling predicts 
no metallic state at all. To get a chance to describe a 
metallic state, one needs to add more ingredients, e.g., 
electron-electron interactions, to the theory. Once we 
deal with a problem which is richer than localization in 
quenched disorder, there are no reasons to assume that 
the one-parameter scaling still holds (see, e.g., Ref. [p0|). 
Therefore, a temperature-independent critical resistivity 
is an assumption rather than a law of nature. 

As an example, consider the following density- and 
temperature-dependences of the resistivity (which do not 
follow from any of the existing theories, but do not con- 
tradict to any of the commonly respected general princi- 
ples either): 



p{n,T) = p (n) + pi(n) 



T 



exp 



Tg(n) 
T 



(1) 



It is natural to define the critical concentration, the 
concentration at which T (n) changes sign: T (n c ) = 0; 
for n > n c (metal), Tb(n) is negative, whereas for n > n c 
(insulator), Tq{ji) is positive. It is also assumed that 
Po(n), pi(n) and Ti(n) are smooth functions of the den- 
sity in the vicinity of the transition point n — n c . Ac- 
cording to Eq.(l), exactly at the critical point (n = n c ) 



T 

p{n c ,T) = p a +pi ( — 



(2) 



As we see, the critical resistivity is temperature- 
dependent and can even diverge as T — > 0. At the same 
time, the resistivity saturates at p = po(n) for T < —To 
in the metallic regime. (In our example, p(T) has a 
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maximum at T = —To /a. We do not think that such 
a maximum is a mandatory feature of the MIT in two 
dimensions.) We conclude that, at least in this example, 
a temperature-independent resistivity is a signature of a 
metal rather than of a critical point. 

Regardless of this particular and rather artificial exam- 
ple, we notice that it is neither easy nor a straightforward 
task to determine the critical point of an MIT. Quan- 
tum phase transitions are zero temperature phenomena, 
whereas experiments are performed at finite T . There- 
fore, to determine the critical concentration, one does 
not have another choice but to analyze the data taken at 
the lowest accessible temperatures. This analysis should 
prove that piT) indeed acquires insulating exponential 
behavior as soon as the concentration gets lower than the 
apparent value of n c . 

We find it more appropriate to use the onset of the 
insulating exponential behavior in p(T) rather than the 
vanishing of the derivative dp/dT for an experimental 
identification of the critical point. We are not aware of 
any reasons to assign the meaning of the critical point 
to a density n* at which the temperature dependence 
of the resistivity is least pronounced, as it was done by 
KK. Indeed, direct measurements of the two quantities, 
n* (from the temperature independent "separatrix" ) and 
n c (from the disappearance of the non-Ohmic hopping 
behavior) confirm their systematic difference for high-/! 
samples: n* is larger than n c by 1-5% (see, e.g., Fig. 3 in 
Ref. [|3C[]). A similar difference arises also from the com- 
parison of critical behaviors in the thermoelectric power 
and conductance [ [3l| . Note that four out of five electron 
densities in Ref. [ |l5|| , n — n\ — 774, fall into this ambigu- 
ous interval. 

III. EXPERIMENTAL DETERMINATION OF 
THE CRITICAL POINT. 

We now turn to the experimental data of Ref. [ [l5"| . 
The authors assume that n c equals to 7.25 x 10 10 cm -2 
(curve 3 in Fig. |l|). According to the definition of the 
critical point, proposed in Sec. [n|, this assumption im- 
plies that at lower densities (n = ni = 7.17 x 10 10 cm -2 
and n — n\ = 6.85 x 10 10 cm~ 2 ), the resistivity increases 
exponentially as temperature decreases. It turns out that 
one can fit neither p(T,n = 71-2) nor p(T,n = n{) with a 
simple exponential dependence. However, p(T,n = n\) 
can be satisfactorily approximated by 



/0 cxp (T, n) = p* m exp [(T 0m (n)/T) r 



(3) 



with either m — 1/2 or m = 1/3 (dashed line A\ and solid 
line B\ in Fig. [l|, respectively). The quality of the fit can 
be seen from the deviations of the experimental data from 
the straight line in the inset to Fig. □]. Although curve 2 
in Fiq. |l| does not seem to behave exponentially, one can 
still try to fit it by function (|^) in two different ways: 

i) one can assume that To m is proportional to \n — n c \. 
(This dependence was observed in previous studies 



of the MIT in Si MOSFETs [ ||-[||). Then, us- 
ing the value of To m for n = n\, one determines 
Tq p for n — 712 by linear extrapolation, whereas p* m 
is found by fitting the experimental data. The re- 
sults of this procedure for m = 1/2 are shown in 
Fig. [j] (dashed-and-dotted line C2). It is quite clear 
that such a fit does not work (the same is true for 
m = 1/3, not shown), 
ii) Alternatively, one can simply make the best fit of 
curve 2 by function (||) for m = 1/2, 1/3, treating 
To m and p* m as independent parameters. Results 
of this fit are depicted in Fig. [j] by curves Ai and 
£>2, respectively. It is customary to assume that 
in the vicinity of the critical point, To m scales as 
\n — n c \ zu . Given the values of To m for n — nx,2, 
one can follow the assumption of KK that n c = 773 
and find 



[lnT 0m (ni)/r 0m (7i2)] _ 
ln[(n 3 - ni)/(n 3 - n 2 )] 

2.23 for m = l/2; 
3.57 for m = 1/3. 



(4) 



It is quite alarming that these values of the criti- 
cal exponent are substantially different from zv = 
1.2 — 1.6 found in previous studies [ 33 ,34]. 

It is also often assumed that in the vicinity of the crit- 
ical point p(T, n) should be symmetric with respect to 
reflection about the separatrix [ §4-36 . To check if this 
relation works, we reflect curve 1 {\8n\ = 713 — 712 = 
0.4 x 10 10 cm~ 2 ) about the would-be separatrix (curve 3). 
The reflected curve is shown by semi-open circles. Had 
reflection symmetry worked, the reflected curve must 
have been located in between curves 4 {Sn — 0.32 cm~ 2 ) 
and 5 (6n = 0.6 cm -2 ), closer to curve 4- I n £&ct, the 
reflected curve crosses both curves 4 and 5. Thus by this 
criterion 77,3 is not a critical density. 

Unfortunately, for n — 712 KK present the data only 
down to 80 mK, whereas for the rest of the densities, both 
smaller (77,1) and larger (773,774,775) than 772, is shown 
down to 35 mK. If observed, a substantial increase in 
p{T, 772) with cooling from 80 mK to 35 mK would pro- 
vide KK with a strong argument in favor of 773 being 
close to the critical density. 

Within the assumption that the data presented in 
Fig. [j] can at all be interpreted within a paradigm of 
a quantum phase transition, the contradictions demon- 
strated above suggest at least one of the following three 
conclusions: 

i) the separatrix corresponds not to density T7 3 but 
to some smaller density; 

ii) the critical exponent is not universal for a given 
material but depends on the sample preparation, 
geometry, etc.; 

iii) the exponential increase of the resistivity at den- 
sity 77,2 is suppressed due to overheating (we will 
discuss this option in Sec. in more details). 
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Alternatively, one may conclude that the data of Ref. 
[ \ldjj do not provide evidence for a quantum phase tran- 
sition between two distinct states of matter, but rather 
indicate a crossover between two different types of (rela- 
tively) high temperature behavior. 

IV. RESISTIVITY OF SI MOSFETS IN THE 
CRITICAL REGION 

As it has already been mentioned, any scenario of a 
quantum phase transition implies that for the densities 
close to the critical point n c the resistivity p(T, n) at suffi- 
ciently high temperatures demonstrates a critical behav- 
ior. This means that the deviation of p(T, n) from the 
separatrix p(T, n c ) is small for small values of the crit- 
ical parameter \n — n c \ zu /T. Keeping this in mind, we 
now recall the main results of earlier experiments on the 
resistivity of high mobility Si MOSFETs in the critical 
region and compare them with KK's results. 

Consider, for example, the resistivity of a sample with 
Mpeak = 55, 000cm 2 /Vs, reproduced from Ref. [ 30 |37j] in 
Fig. |^. Assume for a moment that metallic and insulating 
behaviors can indeed be distinguished by the sign of the 
derivative dp/dT, as suggested by KK. At first glance, 
the application of this criterion to the region T < 1 K 
in Fig. H is rather unambiguous. Indeed, n* determined 
from the condition of T-independent p corresponds to the 
dashed separatrix. 




T (K) 

FIG. 3. p(T) over the range from 0.016 to 3K. Dif- 
ferent symbols correspond to n from 0.449 to 0.989 
(in steps of 0.054), and further 1.1, 1.2, 1.42, 1.64, 



1.74, 2.82, 3.9, 
Refs. [30,371. 



4.98 x 10 cm . Reproduced from 



On the other hand, using the onset of activated and 
non-Ohmic conduction as a criterion for the transition, 
we find that the critical density n c is lower, n « 0.719 x 
10 11 cm -2 and corresponds to the "tilted" (6th) curve 
from the top. At higher temperatures, some of the 
"metallic" curves change their slopes and follow quite 
closely those "insulating" curves, which correspond to 
lower densities (see solid curves in Fig. ||.) The region of 



the {T, n} plane, where such a behavior occurs, is defined 
as the critical region. Certainly, the density range of the 
critical region is narrower for lower temperatures. 

As a rule, the resistivity depends on the temperature 
in this region. This feature, which can be found in al- 
most any published data, is demonstrated in Figs. || and 
|3|. For example, p(T, n) in Fig. ^ increases by approxi- 
mately a factor of 2, as T changes from 3 K to IK, both 
for n — 0.719 x 10 11 cm~ 2 (the 6th curve from the top) 
and n = 0.773 x 10 11 cm -2 (the 7th curve). 

In fact, the p(T)-dependence at higher temperatures 
is usually more pronounced in higher-mobility samples. 
This is illustrated by Fig. ||, which is reproduced from 
Ref. [B3I. 
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FIG. 4. Resistivity vs T near n c for 3 samples with 
mobility decreasing from top to bottom. Reproduced 
from Ref. [32]. 
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The slope of the dashed line for the sample with 
Mpcak = 71,000cm 2 /Vs (Fig. ||a) is negative and its ab- 
solute value is about 10 times bigger than in Fig. [|b 
(Mpeak = 33, 000 cm 2 / Vs). The reduction of the slope 
with /ipeak persists down to relatively low mobilities un- 
til eventually the slope changes sign. This is seen from 
Fig. ||a, reproduced from Ref. [ 37 1, where the authors 
presented p(T) for a sample with /i poa k as small as 5, 000 
cm 2 /Vs. 




0.16 



0.15 



0.14 



FIG. 5. Resistivity vs temperature for a low mobil- 
ity sample (Si-39). Numbers at the curve indicate the 
densities in 10 10 cm -2 . Reproduced from Ref. [37]. 



Even for this low-mobility sample, the /^In- 
dependence at the highest density measured (n = 54.4 x 
10 11 cm -2 ) shows no signs of saturation in the accessible 
temperature region (above 100 mK), see Fig. [s]b. There 
is no consistent theoretical understanding of this /x pea k- 
dependence of the slope yet. Such an understanding will 
probably come together with the description of the "crit- 
ical resistance" , which apparently changes from about 
0.3e 2 /h at low mobilities to several units of e 2 /h at high 
A*pcak[ |3^j3^, (see also discussion in Refs. [ p8| , p9[ ). On 
the other hand, it is natural that for the intermediate 



mobilities, the slope is small, because it should vanish 
somewhere in the range from 5,000 to 30,000 cm 2 /Vs. 
The mobility of KK's sample is 27,000 cm 2 /Vs, i.e., right 
in this range. Therefore one should not be too surprised 
by the fact that in some range of the concentrations the 
temperature dependence of the resistivity is rather flat, as 
it was observed by KK. 

Last but not least: a variety of the p(T)- dependences 
in the critical region indicates by itself the absence of uni- 
versality in observed MITs [ [?0| and raises doubts that Si 
MOSFETs undergo a genuine quantum phase transition. 



V. "HOW COLD ARE THE ELECTRONS?" 

Heating of electrons by the applied source-drain field 
as well as by external noise is a common problem in 
low-tempcrature transport measurements. This problem 
turns out to be especially serious in Si MOSFETS in the 
vicinity of the "metal-insulator transition" and at tem- 
peratures < 100 mK due to *) weak temperature depen- 
dence of the resistivity in this range of densities and **) 
weak electron-phonon coupling. As a result, the inter- 
pretation of resistivity measurements becomes rather am- 
biguous. For example, the temperature interval, in which 
the resistivity appears to be temperature- independent, 
may look to be much broader than it actually is. The 
goal of this section is to demonstrate the seriousness of 
this problem. 

In macroscopic Si MOSFETs, as well as in bulk Si, 
electrons couple to phonons only via the deformation po- 
tential. As this coupling is rather weak, the electron tem- 
perature, T e , may exceed substantially the bath temper- 
ature, Tbath! it is also not easy to control T e . To provide 
convincing data, one has to use a "thermometer" capable 
of measuring T e directly. Any observable which depends 
on T e , other than the (zero-magnetic-field) resistivity it- 
self, can be used for this purpose. In earlier studies on 
high mobility Si-MOSFETs, Refs. [ which 
reported results of the resistivity measurements down to 
Tbath — 18 — 25 mK, the electron temperature was at- 
tempted to be measured via 

i) the amplitude of Shubnikov - de Haas (ShdH) oscilla- 
tions, 

ii) temperature dependence of the hopping resistivity, 

iii) linearity of I — ^-characteristics, and 

iv) T e -dependence of the phase relaxation time. 
Unfortunately, results on T e in the up-to-date measure- 
ments of 2D MIT are neither fully consistent with each 
other nor convincing enough below Tb a th = 250 mK. 

No independent measurement of T e is reported in Ref. 
[ [l^] . It is only mentioned that the source-to-drain bias, 
U, was less than 200 //V "to ensure that the total power 
dissipated in the sample was less than 10 -13 W" . As it is 
not clear whether this is enough to prevent electrons from 
being overheated, it is worth discussing electron heating 
in Si MOSFETs. 
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Electrons are driven out of the equilibrium by the ap- 
plied voltage and/or by external noise. Let P be the 
power deposited in the electron system. In the stationary 
state, all this power leaves the system either with elec- 
trons (through contacts) or with phonons. The phonon 
mechanism dominates at higher temperatures. If T is low 
enough, this mechanism can be neglected compared to 
electron out-diffusion, in which thermal conductivity of 
the electrons determines the heat balance. As our task is 
to estimate heating at T ^lOOmK, we first discuss what 
happens without phonons. 

Given the temperature increment AT = T e — Tbath ; 
the power, which is carried out by electrons through the 
leads, can be estimated as (see Refs. [ Q and also [ Q 
for recent discussion): 



screened deformation potential. The corresponding re- 
laxation time for Tbath *C TiUfS = 3.65-\/re [K] (where s 
is the speed of sound, and n in 10 11 cm -2 ) can be written 
as [ HH 



120 



7"c— ph 



E 2 m*k 2 F a 3 t T bath \ (Th 



h 2 Ms V k FS J 



oath 



(9) 



where S is the deformation potential constant, M is the 
atomic mass, ao is the effective lattice constant defined 
in such a way that M/clq is equal to the mass density, 
and k is the inverse screening length. Comparing (||) to 
(B), one finds that out-diffusion dominates, provided that 
Tbath is smaller than 



P = 



TaathAT 

R ' 



(5) 



where R is the resistance which differs from the resistivity 
by the aspect ratio. To obtain a lower bound estimate 
of the bath temperature Th at which heating becomes 
strong, i.e., AT ~ Th, we neglect heating due to external 
noise and assume that the main reason for heating is the 
source-to-drain bias, U. Accordingly, P = U 2 /R. Ex- 
pressing the ratio AT/Tb a th through P via Eq. (||), we 
obtain 



AT / 
= PR 1 



T 



bath 



\2nTi 



bath 



-(.. 



eU 



V 27rT ba th 



(6) 



Strictly speaking, Eqs.(|5j,|6|) are valid only when AT <C 
Tbath- Nevertheless, one can use them to estimate Th as 



T h [mK] = 1.8U[pV}. 



(7) 



For U ~ 200pV, Eq. (Q) gives T h = 360 mK. For lower 
bias, e.g., for U = 50 pV this estimate gives Ti = 90 mK, 
whereas the experimental dependence T e (P) [ Q taken 
at Tb a th = 140 mK shows that at this bias electrons are 
already overheated: T e = 1.4Tb a th- Despite its simple 
form, Eq. (0) is quite general and universal. For exam- 
ple, it is also valid in the presence of Fermi-liquid inter- 
actions, which do not affect the Wiedemann-Franz law 
[@- 

Taking into account the phonon mechanism of energy 
relaxation does not change our estimate of characteris- 
tic temperature Th substantially. This is so because Th 
is low enough for out-diffusion to dominate over phonon 
emission. Indeed, the electron-phonon energy loss rate 
for weak overheating can be written as 



r> 71-2 T-w^bathAT 

<J 7~c — ph 



(8) 



where vp is the density of states at the Fermi level, L 
and W are the sample length and width, respectively, 
and r e _ph is the electron-phonon relaxation time. We 
assume that electrons are coupled to phonons via the 



-hk F s—£Vr 2 
i pe* 



£ F 



(10) 



where £ = Ms 2 (fi 2 /m*al) /S 3 and V = h/m*sa . Sub- 
stituting material constants for a Si MOSFET into (pT 
one obtains 



T*[K] = 0.3 (pn 2 y 1/5 L- 2 ^ 



(11) 



where p, n, and L are measured in h/e 2 , 10 11 cm -2 , and 
mm, respectively. Using experimental results on electron 
heating in Si MOSFETs [Q, one estimates the observed 
crossover temperature as 0.7 K, whereas Eq. ( pT| ) gives 
T* = 0.6 K for the experimental conditions of Ref. [ 45|. 
For the conditions of another experiment on heating [ 49 1, 
Eq. ( pl| ) gives T* — 0.3 K, which is just the lowest tem- 
perature of this measurement. No clear crossover from 
the phonon to out-diffusion mechanisms was observed in 
Rcf. [ H for T > 0.3 K (although the T-dependence of 
the energy-loss rate does start to slow down at T ~ IK). 
This is again consistent with our estimate. We thus con- 
clude that Eq. (|ll|) is consistent with the experimental 
data and can serve as at least a lower bound for T*. 

For typical values of p ~ e 2 /h, n ~ 10 n cm~ 2 , and 
L = 0.1 - 1 mm, Eq. (fnj) gives T* ~ 0.3 - 0.75 K, which 
is of the same order as Ti estimated above as 0.36 K. 
Note that in deriving ( pT| ) we assumed that p was T- 
independent. Taking the metallic-like p(T)-dependence 
into account enhances out-diffusion at low temperatures 
and thus shifts T* towards even higher values. Also, 
taking into account possible external noise would further 
increase the value of Th. The lowest bath temperature in 
Ref. [ |T| is T b = 35 mK. Using Eq. (Q), we see that in 
order to prevent electron heating, the bias voltage has to 
be much smaller than 20 pV, which is ten times smaller 
than the upper boundary for U indicated in Ref. [ ^] . 

We thus see that there are serious reasons to doubt 
that the electron temperature in KK's measurements was 
below ~ 100 mK. 
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VI. TEMPERATURE-INDEPENDENT 
RESISTIVITY 

Now it is time to analyze the first of the two main ar- 
guments that KK brought in favor of the "true quantum 
phase transition" , namely, the temperature-independent 
resistivity at what they believe to be the critical point. 
More specifically, KK claim that they observed no 
T-dependence (within a 10% margin) in the interval 
35 mK < T < IK. They analyzed the data in terms of our 
recent theory of Anderson localization by temperature- 
dependent disorder [ || and came to the conclusion that 
within this theory a nearly constant p(T) would imply a 
"remarkable coincidence" . 

We begin our discussion of this issue with summarizing 
briefly the argument of KK. Afterwards, we present our 
theoretical counterarguments. 

The theory of Ref. [ ||] describes the T-dependence of 
the observable resistivity, p, in the presence of two fac- 
tors: i) the T-dependence of the Drude resistivity p d , 
and ii) Anderson localization. Each of these two factors 
brings a T-dependence of its own and the resultant p(T)- 
dependence is described by the following scaling equation 



1 



■ 7 



dhxp p 
rflnT = 2 



1 



P(pd 



din 
dlnT ' 



(12) 



where p is measured in units of h/e 2 , /3(x) is the Gell- 
Mann-Low function, whereas p and 7 parameterize the 
phase-breaking time as 



ji-p 1-27 



(13) 



The T-dependence of pa, entering the RHS of Eq. (Il2| ), 
leads to a variety of behaviors in p(T). In Ref. [ |we 
discussed in particular the situation when p{T) exhibits 
a maximum at a rather high temperature T max (close to 
the starting temperature of the RG flow). 

Is a very slow variation of p(T) over a wide T-interval 
permitted in this model? KK answer this question nega- 
tively. Their argument goes as follows: to have dp/dT = 
within some interval of temperature, one has to require 
that the RHS of Eq. ( |l2| ) is equal to zero within this in- 
terval. This is possible only if the Drude conductivity pd 
is related in a specific way both to the /3-function and to 
t v . Such an exact relation is unlikely, given the different 
origin of these quantities. 

This argument sounds to be formally correct. Nev- 
ertheless, as we will demonstrate shortly, it is possible 
to achieve a very small variation of p(T) within a wide 
temperature interval without imposing any constraints 
of this kind on pd- It turns out that the variation of 
p(T) is small enough and the T-interval is wide enough 
to agree with the experiment. This possibility results 
from the fact that there is a whole family of the Pd{T)- 
functions, which are parameterized by the electron den- 
sity and some other parameters, such as the peak mo- 
bility. It is the freedom in fine-tuning these parameters 



that allows one to suppress the variation of the observ- 
able resistivity dramatically. 

For the purposes of this section, we adopt the following 
working definitions, consistent with those used in Ref. [ 
p|: by "metal" ("insulator") we understand the region 
in which dp/dT > 0(< 0) or dp d /dT > 0; by "critical 
region" we understand the region in which p ~ 1. 

As was discussed in Ref. [ |J, a maximum in p(T) 
results from a competition between metallic- like Pd{T) 
and localization effects, the latter being controlled both 
by pd(T) and t v . At higher T, when pd is sufficiently 
large, localization can already be strong enough to en- 
sure the negative sign of dp/dT despite pa(T) decreas- 
ing with temperature. However, as T goes down, pd de- 
creases, localization weakens and dp/dT changes its sign 
from negative to positive, thus a maximum in p(T) oc- 
curs. This can happen provided that certain conditions 
are satisfied. In particular, one has to require that 



Pd(T ) > p c 



(14) 



where Tq is the temperature at which t v is comparable 
to the transport mean free time (the flow of Eq. jl^ ) 
starts at T = Td) and p c is some critical resistivity de- 
termined by a particular form of pd(T). For example, if 
Pd oc T q , critical resistivity p c is a certain function ofp/q 
(see Eq. (9) of Ref. [ §). In this case, there is another 
condition for the maximum in p(T) to occur, namely, 
P > 2q. 

As the resistivity pd(To) depends on the electron den- 
sity n, it can be tuned by the gate voltage. This tuning 
can drive the system between the domains of parameters 
corresponding either to a maximum or no maximum in 
p(T). Let us start at Pd{Ta) > p c , so that the maximum 
in p is at T = T max . As pd approaches p c from above, 
Tnax remains finite, while the maximum flattens out and 
disappears as soon as Pd(T)) approaches p c . At this mo- 
ment, both dp/dT and d 2 p/dT 2 vanish, and the variation 
of p around T max takes place solely due to higher [n > 3) 
terms in the Taylor expansion of p. 

How large is the temperature interval, in which this 
variation does not exceed some given value? Denoting, 

t = InT, 

^rnax — ln^max; 
/?d = /3[Pd(*max)], 

y(t) = lnp d (i), 
P'd = d(3(p)/d\np 

p=Pd 

(% = d 2 l3(p)/d\np 2 
we obtain from Eq. ( |l2] ) 
= 



dp 
~dT 
d 2 p 
dT 2 



T=T TI , 



T=T m 



. , , \ 'P_ Pd 



2 (3d 



(15a) 



= 



P'd 



y{^) = — " d [yfenax)]'- (15b) 

Pd(pd - 1) 



9 



The third (logarithmic) derivative of p is then found from 
Eq. (O) to be 



d 3 In/ 



dt 3 
where 

and 



s ^ 2/? d 2 - p d p> d > 



= {II A) 



Pd - 1 A 



p d dt 3 



A = l/(3(p) - 7 



Pd 



8(/3 d - I) 4 



Pd 

{2p' d 2 -p d p';yp d -\ 



c)j/(*m 



(16) 

(17) 

t) (18a) 
. (18b) 



(In going from (18a) to ( |18b| ), we used (15a, 15b)). 

For a sufficiently narrow interval of t around t max , it 
suffices to keep only the lowest (cubic) term in the Taylor 
expansion of lnp(T) as a function of ln(T/T max ): 



In 



Pit) 

p(^max) 



= (1/6A) 



S- 



ffd-l d 3 lnp d | 
/3 d dt a !*=*» 



^ (^ ^max) 



(19) 



As y is supposed to be a smooth function of t, one can 
estimate d n y/dt n as y/t™, where y is the typical value of 
y in the interval i. In the "critical region", p ~ 1, which 
means that y ~ 1 as well, and hence /3d — ^ Pa — 1- 
Thus the two terms in the square brackets in (|19| ) are of 
the same order and (19) reduces to 



lnp(i) - lnp(i n 



(S/6A) (t - t n 



(20) 



Let us now estimate the number D of decades in temper- 
ature 



Thigh /Tiow = 10 D , 



for which 



In /»(*)- La p(W) I < 10%. 



(21) 



(22) 



(KK found that their p is T-independent with this 
accuracy.) Using an interpolation formula P{p) — 
— In (1 + 2p/n) [ and choosing p d (T max ) = tt/2, we 
have 



D <2.S\A\^ 3 /p. 



(23) 



For p = 1 and |A| = 1, we get D < 2.3, i.e., the resistiv- 
ity stays within a 10%-margin for more than two decades 
in T. (Note that a numerical coefficient of the order of 
unity, which we neglected in deriving (ppf), would enter 
( p3| ) only under the cubic root, and is thus unimportant.) 
Fig. H demonstrates p(T) calculated from Eq. ( |l2] ) for p d 
taken from the model of Ref . [ ^l| : 



Pd = Pi + p(T/T ) q 

(1 + \6\T /cT) q , for S < 0, "insulator"; 
e -5T /T ; for g > Qj « meta i" ; 



(24) 



where J = <5o(™ — n c )/n c characterizes the distance from 
the critical point and c ~ 1. The dimensionless coeffi- 
cient Sn depends on details of the model (cf. Eq.(ll) of 
Ref. JJ5|). For the purposes of the present paper, we view 
Eq. Q24 ) as simply a phcnomenological form for p d (T), re- 
gardless of the model [ ^l[ it was originally derived from. 
This form is consistent with the experiment, if 5q — 1. 

The value of S — corresponds to the central line in 
Fig. [| which on this scale is essentially constant for more 
than two decades. Zooming in (cf. the inset in Fig. @), 
one sees that in fact p(5 — 0, T) stays within a 10%- 
margin over exactly three decades in T . 




FIG. 6. Solution of Eq. (|T|) for different val- 
ues of 5 (distance from the "critical point"). 
p = l,q = 1/4,7 = l.Pi = 0.28, p = 2.3. Inset: 
the central (almost flat) curve (8 = 0) over a wider 
temperature range. 
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KK claim that the observed p(T) stays within this mar- 
gin over only 1.5 decades, which is smaller than even 
the conservative estimate (|2^). (As it was discussed in 
Sec. 0, a more realistic estimate of KK's T-interval is 
one decade). A decrease in S by 1% (S = —0.01) leads 
to a weakly insulating behavior: p(T) increases by 30% 
as T decreases by 1.5 decades. An increase in S by 4% 
(6 = 0.04) leads to a well-pronounced metallic behavior: 
p(T) drops by about a factor of two over 1.5 decades in 
T. These features are in quantitative agreement with the 
experiment [ |l5| . (For the sake of simplicity we assumed 
that parameters p\ and p do not depend on the density, 
i.e., on S. Taking these dependences into account should 
lead to even better agreement with the experiment.) 

We emphasize that no fine tuning between the /?- 
function and pd was used. If it is the competition be- 
tween the decrease in Drude resistivity with temperature 
and the localization effects that leads to the maximum in 
the observable resistivity, then the /9(T)-dependence flat- 
tens out, as the density approaches the threshold for a 
maximum. It turns out that if the exponent p is not too 
big (which is usually the case) , this flattening suppresses 
the variation of p(T) below the experimentally observable 
level over several decades in temperature. Even an over- 
simplified model [ H which does not take into account, 
e.g., electron-electron interactions (except for as a pos- 
sible phase-breaking mechanism), can easily produce al- 
most constant resistivity in a temperature interval, which 
is two orders of magnitude wider than the experimental 
one. Therefore, we do not find it too much of a "re- 
markable coincidence " that at some density the system 
demonstrates a fairly constant p. 

We now turn to the data presented in Ref. [ [L5| . 
The authors emphasize strongly that at n = = 
7.25 x 10 11 cm -2 the resistivity is almost temperature 
independent in the temperature range 35 mK - IK. They 
also add that p{T) decreases with T for T > 1.8 K. It 
would be interesting to know what happens at interme- 
diate temperatures: IK < T < 1.8 K. This point being 
unclear, KK fill in the missing temperature interval by 
combining their result with the one obtained on another 
system. Indeed, they write: 

"In combination with the results of Ref. [10] (Ref. [ |5^] 
of this paper- AMP) where the temperature-independent 
curve (with essentially the same value of resistivity) was 
observed in the temperature range 250 mk-1.8K in an- 
other 2D system in silicon, we further allege that there 
is no observable T-dependence at n s — n c in the tem- 
perature range 35mk-1.8 K...". (Ref. [ [[5), p. 3, second 
paragraph.) 

Parenthetically, we note that the resistivities of combined 
curves differ by 35%. As far as we understand, tempera- 
ture intervals of experiments on different systems are not 
additive parameters, and thus the procedure described 
above is not justified. 

Summarizing this part of our discussion, we take the 



liberty to describe the experimental situation in the fol- 
lowing way: 

In a sample from the intermediate peak mobility range 
(where temperature dependence of p is known to be 
quite weak at high T), one can find a density such 
that p does not change for more than 10% within about 
one order of magnitude in T. This statement is a re- 
sult of measurements on a single sample in the interval 
100 mK< T < IK (when a realistic estimation of the 
electron heating is taken into account). (It will be two 
samples if results of Ref. [ |52| for 250mK< T <1.8K are 
taken into consideration.) 

In our opinion, there are neither theoretical nor exper- 
imental reasons to believe that the density corresponding 
to the weakest p(T)-dependence is indeed the critical one. 
i.e., that p(T — > 0) — > oo and p(T 0) is finite for lower 
and higher densities, respectively. The lack of a pro- 
nounced temperature dependence of p in a limited range 
of T does not signal any remarkable phenomenon and is 
well-described by a simple theory of Anderson localiza- 
tion in a temperature-dependent disorder (Ref. [ [5) ) . 



VII. ON THE APPARENT LINEAR 
TEMPERATURE DEPENDENCE OF p(T) 

We now turn to point (2) of Ref. [ |l5| )] regarding the 
functional form of the observed p(T)-dependence. Ac- 
cording to KK, this form is a) linear and b) qualitatively 
distinct from those observed by other authors. Fig. 
reproduces Fig. 3b of Ref. [ , where the "almost lin- 
ear dependence" is demonstrated. The data in Fig. fj] is 
data 5 of Fig. |l|, displayed in truncated (T < 400 K) and 
full (T < 1200 mK) intervals. The truncated interval, 
in which p(T) is supposed to be "almost linear" , corre- 
sponds to Fig. 3b of Ref. [[f§. 

We begin with an attempt to fit untruncated data 5 of 
Figs, [j], [7] in the whole interval T < 1.2 K by an empirical 
expression [ Q 

p(T) = p +p 1 exp(-'^j , (25) 

which is a simplified version of Eq. (1) with a = 0. 
This expression is known to work reasonably well for 
the metallic phase not only in Si MOSFETs [ ||J5|] but 
also in other 2D systems exhibiting the MIT-like behav- 
ior[||. 

As is seen from Fig. 0, this attempt is quite success- 
ful. We achieved more than just acceptable fit in a much 
broader temperature interval than the one in which p(T) 
looks "almost linear" . Of course, it is not a much of an 
achievement to fit a smooth curve by a function with 
three free parameters (to demonstrate a flexibility in 
choosing the fitting parameters, Fig. [?] presents two sets 
of parameters that can be used). On the other hand, we 
see no reason to argue that the p(T)-dependence of Ref. 
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[ [15| at n — ri4 and n — n 5 (see Fig. 1) is qualitatively 
different from those observed in other experiments. 
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FIG. 7. Temperature dependence reproduced from 
Figs. 1, 3 of Ref. [15]. Region to the left of the verti- 
cal dashed line is the truncated T-range in Fig. 3a of 
Ref. [15], in which p(T) is claimed to be linear. Bold 
dashed and continuous curves: fits of the data with 
Eq. ©. 

Moreover, such "almost linear" /o(T)-dependences are 
typical for low-mobility samples in the "metallic" regime. 
For instance, compare Fig. [7] with Fig. |5|b, in which p(T) 
for a low-mobility sample is plotted over a much wider 
temperature range (90 mK - 4K), and for n about 50% 
higher than what can be called a critical concentration. 
The similarity of these two plots is quite clear. 

We find it quite unlikely that the quasi-linear p(T)- 
dependences of Figs. [l] and fj] (as well as Fig. 5) manifest 
an exotic non-fermi-liquid state. Taking into account the 
whole set of features of the metallic state in Si-MOSFETs 
(which we listed in Ref. [§), we conclude that these de- 
pendences must have a classical explanation (see the next 
paragraph for our definition of this term). For example, 
one cannot exclude that such behavior is caused by the 
temperature dependence of the screening length 
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or by the temperature dependence of the charge traps' 
population in the oxide [ ^l| . 

We conclude this discussion by remarking that there 
seems to be some confusion in recent literature on MIT 
in 2D (see, e.g., Ref. [ ^6|) regarding the distinction 
between "quantum" and "classical" behaviors. We be- 
lieve that this confusion is mainly semantic. We call a 
regime "classical" as long as electron transport can be de- 
scribed in the framework of the Boltzmann theory. This 
does not necessarily mean that electrons obey Boltzmann 
(non-degenerate) statistics -in fact, the Fermi-liquid the- 
ory was originally formulated in terms of the Boltzmann 
equation [ p9|. It also does not exclude that a scattering 



cross-section is described in terms of quantum mechanics. 
(A more common name for this regime is "semi-classical" 
[ Q , although the semi-classical approximation may not 
necessarily be employed for the calculation of a scatter- 
ing cross-section.) However, quantum-mechanical effects 
arising due to interference of electron waves, e.g., weak 
localization, cannot be described by the Boltzmann equa- 
tion. We call the regime "quantum" , if these effects are 
important. Although we agree with the authors of Ref. 
[ [16| in that "the formation of a Fermi surface" is an "in- 
trinsic quantum-mechanical effect" , this does not imply 
that the metallic-like temperature dependence of the re- 
sistivity, whose onset is correlated with the Fermi-surfacc 
formation, is of the quantum nature, in a sense of the def- 
inition given above. 



VIII. AT WHAT TEMPERATURES WEAK 
LOCALIZATION MAY BE OBSERVED 



There is one more point to discuss regarding the metal- 
lic behavior of Si MOSFETs. In addition to the almost 
linear dependence with dp/dT > 0, KK emphasize that 
no signs of weak localization (WL) was observed, al- 
though they went to temperatures as low as T ~ 0.008Tf, 
where Tp is the Fermi temperature. Giving the authors 
the benefit of the doubt, we assume the most favorable 
scenario in which electrons are not overheated. Neverthe- 
less, under this assumption the apparent absence of WL 
is not as surprising as KK present it. In fact, we believe 
that even the bath temperature was too high to observe 
the localization upturn in the T-dependence of the resis- 
tivity. Indeed, experimental studies of high-mobility Si 
MOSFETs in the wide range of electron densities [ p2 61 
show that WL effects are well-pronounced in the (perpen- 
dicular) magnetic field dependence of the resistivity. The 
observed magnetoresistance is in a quantitative agree- 
ment with the conventional theory [ , even for densi- 
ties rather close to the critical one. On the other hand, 
even for n> n c observation of the WL effect in the zero- 
field p(T)-dependence [ f42[,[s| (as well as the effects of 
electron-electron interactions) requires a precision of the 
order of Sp{T)/p(T) < (l/n){e 2 /h)p ~ 0.1%. 

It turns out that the observed dependence in the range 
of "high temperatures", T = (0.1 — 0.5)Ep, is too strong 
to be explained by WL and quantum interaction effects, 
and can be attributed to the T-dependence of the Drude 
resistivity. However, the slope \dp/dT\ decreases with 
T. For relatively large densities, one can reach the tem- 
perature region, T < Tq, in which the slope becomes 
comparable to that of the WL correction, \dp / dT\w l ~ 
(e 2 /ir)p 2 T. In this region, quantum interference should 
contribute significantly to the p(T)-dependence. Authors 
of Ref. [ |qj42| found an empirical relation: Tq < 0.007 Ep 
(see Fig.7^))T As Ep is proportional to n, it becomes pro- 
gressively more challenging to cool electrons below Tq as 
n c is approached. Ep ~ 5.5 K for the highest electron 
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density reported by KK and thus Tq = 40 mK. There- 
fore, in order just to reach the high temperature edge of 
the quantum transport region, where WL effects are seen, 
electrons should be cooled below 40 mK. The localization 
upturn may be expected to occur only at temperatures 
substantially smaller than even this one. In other words, 
the lowest temperature reported by KK (35 mK) is still 
too high for localization effects to be observed. 

Why WL manifests itself so differently in the mag- 
netoresistance and in the p(T)7 This can be explained 
naturally by assuming that the Drude resistivity has a 
pronounced T-dependence, which masks quantum cor- 
rections. On the other hand, pd is not expected to vary 
substantially in weak magnetic fields, thus WL can be 
seen in magnetoresistance. 



IX. CONCLUSIONS 

Some concluding remarks: 

i) the apparent metal-insulator transition in 2D is a 
very interesting and unexpected phenomenon. Al- 
though we have limited our discussion to Si MOS- 
FETs, which exhibit the strongest, among other 
2D systems, "anomalous" metallic behavior, there 
is a whole variety of interesting effects not only in 
transport but also thermodynamic properties, e.g., 
compressibility [ ^3j,|i)), observed in various 2D 
heterostructures. Whether all these observations 
have a universal explanation or there is a number 
of different, system-specific mechanisms at work, 
remains to be seen. 

ii) we do not believe that the experimental results 
accumulated up-to-date provide a convincing evi- 
dence for this phenomenon being a true quantum 
phase transition between two distinct states of mat- 
ter. More precisely, there is no evidence that the 
"metallic phase" (as defined by the sign of dp/dT) 
is a new state of matter. On the contrary, a whole 
set of features-see Refs.[ ||,^|-is consistent with a 
conventional behavior of a disordered Fermi liquid. 
This statement does not negate i). There is still 
no consensus on the origin of the effect, and this 
problem requires the most serious and intensive in- 
vestigations. 

iii) For a successful resolution of the problem, it is 
probably not sufficient to concentrate on the tem- 
perature dependence of the resistivity in the vicin- 
ity of the transition. Such features as a weak tem- 
perature dependence of the resistivity at a certain 
electron density in a limited (though large) tem- 
perature interval allow for different interpretations 
and, unfortunately, can not provide an unambigu- 
ous information about the zero-T state of the sys- 
tem. In particular, it is impossible to determine 



even the critical density without making rather ar- 
bitrary assumptions. 

iv) There is a serious experimental difficulty that up to 
now prevented a substantial increase of the temper- 
ature interval, in which the resistivity can reliably 
be measured in Si-MOSFETs. It turns out that in 
existing samples phonons can hardly cool the elec- 
tron gas below a fraction of a Kelvin. On the other 
hand, the systems are quite noisy, and thus the ap- 
plied voltage cannot be reduced much below the 
currently used level. As a result, even in the ab- 
sence of non-equilibrium noise, one cannot neglect 
electron heating at T ~ 100 mK. 
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